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Transient Modeling of Hybrid Rocket Low Frequency Instabilities 

M. Arif Kjnbeyogli*, SfcneDe Z3wa f , Brian Cantwell 1 and Greg Ziffiac* 
Abstract 

A comprehensive dyaanic mode] of a hybrid rocket has been deve l oped io order to understand and 
predict the transient behavior including instabilities. A linearized version of the transient model predicted die 
low-frequency cham b e r p re snuc oscillations that arc commonly observed m hybrids. The saner of the 
instabilities is based on a complex coupling of thermal transients in die solid fuel, wall heat transfer blocking 
due to fod i c gcaio n rate and the tra a si c jits ia the bonuday layer that forms on the fuel surface The 
oscillation frequencies predicted by the linearized theory arc in very good agreement with 43 motor test 
results obtained from the hybrid propnhicn fterte. The motor test resate used a die co m p arison cover a 
very wide spectrum of parameters including: 1) four separate research and development programs, 2) three 
diffe re nt otitfizm (LOX, GOX. N*0), 3) a wide range of motor dimensions (Leu from 5 inch diameter to 72 
inch diameter) and operating conditions and 4) several fuel formulations. A simple universal scaling formula 
for the frequency of the primary oscillation mode is suggested. 
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V V ' 
r » • 

Port and motor volumes 

z : 

Axial distance along the port 

5. 

Boundary layer thickness 

A P^: 

Chamber pressure oscillation 

r • 

Ratio of specific heals 

k : 

Thermal diffisivny of fad 

Pf* P - 

Fuel and average gas densities 


Bombay layer response time and 

T ' 

thermal lag time 

*#* 

Gas phase response coefficients 

P’ 

Subscripts: 

Average gas viscosity 

a : 

Ambient 

1: 

Fkn pqtnrinrion variable 

e : 

Free stream vahie 

r«f : 

Reference qualities 

s : 

Surface 





The allowable thrust oac ib b o g for an 
o p erati ona l propateoa lyrtem are fau l te d by certa in 
p ractical considerations iodafag aatloai oo 
loads on tbe vriidr stwetwo and abo oe the 
payload. This min wdndin a stability a critical 
issue daring the d evel opment of a near p rop abk i 
system. The combustion instabilities of solid and 
liquid systems have been stndied in depth and 
partially s u ccess ful theories along with certain 
practical rules for designing stable systems hare 
been developed 1 ' 23 . Even with this extensive 
knowledge base, development of stable rocket 

resources firing the d esign nd testing phases of a 
new system. Even though he ti a niicn t op erati on 
and mstabiity mechanisms of hybrids have not 
been e xp l o red as exten sivel y as Ac more mature 
chemical systems, past e xp eri en ce shows that 
hybrids do not have the catastrophic mstaNifirs 
that liquid and solid systems may potentially 
present. 


hi fact, hybrid systems typically show 
finite n y tito d e (Le. 2-20% ims of mean chamber 
pressure) low-freq u enc y cham ber pressure 
osenbrions* 5 * 2 *. This most common hybrid 
instability is in the form of limit cycle oscillations 
with frequencies much than the first 


knpbdhal acoustic mode of die combustion 
cham ber. 


Even though the exact cause of the low 
frequency hybrid instabilities is not yd known, a 
few pla u s ible theories c3Ust 5 ' fJ * JU2 . Most of Ac 
suggested theories lack Ac maAegnric a f f or malism 
that exists in the solid ad Squid instability models. 
The fun da mental difficulty h producing a 
mechanism that would generate instabilities comes 
from tbe fact dot hybrid bum rales are typically 
indqmkn of the rhihrr p re ss u re . This* for 


the thermal lags h the solid and the chamb er gas 
dynamics which is the sourer of the L* i nstab il ities 
in sohd rockets 1 . 


One particular theory, tbe TC coupled 
model 1 ***, which is based on be coupling of 
thermal lags in the solid and the lags in the gas 
phase transport, has predicted (based on a 
math ema ti ca l formalism) fed ana p roduc ti on 
oscillations that would drive the o bs e rve d press u re 
oscillations. The purpose of this paper is to extend 
the TC coupled theory to include the effects of the 
ch a m ber gas dynamics and c o mp a re the theory 
predic t ions with the mchor test data. The merits and 
shortcomings of theTC coupled model wS abo be 
rfisemsed. 

M o la r s; 


Our ap proach to modeling the dynamic 
behavior of hybrids is to isolate Ac subsystems of 
the motor and consider every single subsystem 
individually. For a fuB de scri p ti on of motor 
transients one has to consider the following 
subsystems. 

n) Feed System: In a hybrid motor, the Squid (or 
gaseous) oxidizer needs to be fed in to the 
c omb ustion ih—hrr th ro u gh a feed system. In 
reality the feed system response time is finite due to 
Ac ca pa ci ty of the demons in be system. Accurate 
response de pe nds on be details of the system, 
which is HuJy to be s i gnifica ntly d i fferent for 
every design. For Ass reason we wffl bypass the 
feed system dynamics m our kvestigations by 
assuming tbe caufcg mass flow rate as Ac input 
parameter. Thb is a reasonable assumption for most 
systems bec a use isolati n g element* (Le. ca v il Sing 
v e n t ure, sonic orifice) are present. 

b) Vaporization of the lipte OnAsr 

Complete vap oriz a ti on of Ae oxidizer d ro pl e t s in 
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the combustion chamber requires a certain 
daracterisdc tone dqxntfing cm parameters such 
as the droplet size and the thermal/ flow 
environment of the space where droplet 
vaporization takes place. For the purposes of this 
paper we wfll assume that Ac oxidizer is in gaseous 
phase when it enters the port 

c) Dtffukx and Combustion in the Boundary 
L ayer It re q uires some time for Ae hybrid 
boundary layer properties to adjust to the changes 
in the port velocity or the fbd blowing from the 
surface. The associated dynamics has been 
dnenreed in detail in Reft. 13, 14 aad will be 
briefly covered in section 5 of Ais paper. 

d) Thermal Response in Ae Solid Grain: A 
change in the waQ heat flu to the hybrid fad gram 
can not be followed immediately by fuel production 
due to Ac fantc heat capacity of Ac solid fueL 

e) Chiwfirr Gas Py inirr The chamber pressure 
responds to the changes in Ac mass flow with a 
time scale proportional to the filling time of the 
c ha mber. Acoustic response a also resolved m this 
sub-system. 

h the atenpat sectio ns Ac bet three of 
these snb-systems (seasons c-c) wfll be modele d 
and investigated ridrvidnafly. Enema aBy, Aese 
subsystems wiO be co u p led to give the overall 
system response. The order of mqgniterir es timat e s 
of the time scales of some of Ae important 
processes encountered m a typical hybrid motor ae 
listed in Table 1. In oar modeling, solid aad 
ga se ou s kinetic times are assumed to be fast 
c o mpare d to the other rel e va nt time scales. 


4) Thermal Lag Model: 

The regression rate of the hybrid fuel grain 
cannot respond to Ac changes in the stniacc heal 
flux instantaneously due to Ae finite Aexmal 
conductivity of Ae solid fbd. A complete transient 
model for a hybrid rocket system requires a 
dynamic model to predict Ae regression rale history 
for a given wall heat flux schedule. In order to serve 
this pvpose a thermal lag model that models the 
nonlinear heat conduction m the fuel grain with a 
moving boundary has been formulated and reported 
in Refs. 13 and 14. The schematic of this particular 
model is shown in Figaro 1. 

In Ae thermal bg model, Ae gasification 
and pyrolysis reactions at the surface are boA 
modeled by an exponential expres s ion of the 
Arrhenius type. For chemical reactions, Ae 
exponential constant is an activation energy, 
whereas in vaporization it is the latent heaL In order 
to describe this behavior, wc assign an average 
effective activation energy, E a , resulting in: 

r=A (l) 

Herr T % is Ae s urf ace tem p erature, r is 
the regression rale and it, is Ae lauversal gas 
constant. The mih n mk al founhun and various 
solution techni ques for Ae thermal lag modd are 
presented in Ref 13. One solution of interest is 
obtained by perturbing Ae full nonlinear system 
around the nominal op erating point. The linear 
initial- boundary value problem defined for the first 
order pertwbation quantities has been solved with 
use of the Laplace transformation technique. The 
transfer function between Ae regression rale 
perturbation and Ac ap p l ied heat flax p e rtu rb ati on 
can he written as 


Fr 




2E^j 


fi<*> (» W4*+iJ fTi +2e,e^. 


( 2 ) 


Note that R|(x) and 0,(s) are Ae 
Laplace transforms of Ae wmAraensaond 
regre ssi o n rale pc rtm ba ti ou and the wall heat fht 
pejtubntion, respectively. Here Ae following non- 
dimensional parameters are defined for 
con veni ence. 


* KV.Kf CAT 


*r=(r,Ur m 
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The reference-stale indicated by subscript “reP’ 
corresponds to the nominal ope rati ns point around 
which the system is perturbed. 

Equation 2 relates the wall heat flux to the 
regression rate in the Laplace space and establishes 
one block in our overall transient model that is 
shown schematically in Figure 2. Various 
important characteristics of a linear system 
including stability can be interred from the 
denominator of its transfer function commonly 
referred as the characteristic equation. This 
transfer function for the thermal system, which 
contains a square root term, produces a phase lead 
between the heat flux and the regression rate in the 
low end of the frequency domain 13 . As it will be 
discussed in the following sections, this phase lead 
capability will play a critical role in the production 
of low frequency instabilities. 

5) Gas Phase Combustion Model: 

In the thermal lag mode] we have treated 
the surface heat flux as a parameter that can 
directly be controlled. However, in a hybrid motor, 
the oxidizer mass flow rate is the primary input 
variable that can be controlled directly. The actual 
response of the motor to a change in the oxidizer 
mass flow is rather complicated during a transient. 
As the oxidizer mass flow rate of the motor 
changes, the mass flux at a characteristic point in 
the port reacts to the change as does the turbulent 
boundaty layer developed over the fuel surface. In 
this section we will summarize the modeling of the 
boundary layer combustion dynamics and 
investigate its interactions with the thermal lags in 
the solid. For the sake of simplicity, wc will ignore 
the radiative heat transfer to the fuel surface, which 
is typically a relatively smal l fraction of the total 
heat flux. Unless the radiation dominates the 
convection component of the surface heat transfer, 
the transient model developed in this paper is 
expected to be valid. 

In our preliminary model we assume that 
the boundary-layer response is quasi-steady, 
namely the bound ary-layer diffusion lag times are 
small compared to the thermal lag times in the 
solid. Under the quasi-steady assumption one can 
use the classical approach 15 to calculate the 
response of the wall heat flux to the changes in 
mass flux. However, the results of classical hybrid 
theory cannot be used directly during a transient in 
the solid, due to the fact that the convective heat 
transfer to the wall depends explicitly on the 
instantaneous regression rate through the blocking 
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generated by the blowing of the gaseous fuel from 
the surface. In the presence of the thermal lags, the 
blocking generates a coupling mechanism in the 
gas phase between the regression rate and the heat 
flow to the surface. Based on this understanding, 
the classical theory can be modified to obtain a 
functional relation between wall heat flux and 
oxidizer mass flux-regression rate combination in 
terms of non-dimensional parameters: 

Q c (0 = E l G”K ] - k )R- k ^- k ) (3) 

Note that G 0 = G 0 j{f} 0 ) nf , Q c =Q c jQ r(f and n 

is the oxidizer flux exponent, k is the blowing 
correction exponent first defined by Marxman (i.e. 
C H jC Ho = B k )- Here the local total mass flux that 

appears in the original Marxman formulation is 
replaced by the oxidizer mass flux for convenience. 
The justification for that transformation is given in 
Ref. 14, 

Thus far in the development of the 
transient hybrid combustion theory we assumed 
that the boundary layer responds rapidly to the 
changes in the mass flux compared to die other 
transient time scales in the rocket motor such as the 
gas dynamic lags or the thermal lags. This 
assumption fails to be valid especially for large 
hybrid motors* In order to develop a realistic model 
for the dynamics of the hybrid motor, the boundary 
layer lags must be considered. Since this complex 
dynamic phenomenon is extremely difficult to 
investigate both theoretically and experimentally, 
we consider the simpler cases reported in the 
literature of an incompressible Turbulent boundary 
layer with no blowing or chemical reactions 16 . The 
most important conclusion for this simplified case 
is that the time required for transition from the 
initial equilibrium profile to the final equilibrium 
profile at any axial location, z, is proportional to 
the time of flight of a fluid particle from the 
leading edge of the boundary layer to the specific 
axial location at the speed of the free stream 
flow This very important result can be 
formulated for hybrid boundary layers as 

T bl =c'—. (4) 

Here c is a constant that needs to be 
determined empirically. We will call this time 
required for equilibration, the characteristic 
response time of the boundary layer, r w . It is 
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important to note that the physical nature of the 
boundary layer transient time is not related to the 
propagation of the disturbances with the speed of 
the port velocity as suggested by Equation 4. The 
delay rather depends on the diffusion time scale 
across the boundary layer which is proportional to 
the ratio of the local boundary layer thickness to 

the diffusion speed, r bf = Sji/ m , The diffusion 
speed is defined in terms of the shear stress and 
mean gas density as U* = -fijp - The boundary 

layer delay time, after the substitution of the 
Standard (incompressible) turbulent boundary layer 
expressions 17 for the shear stress, 
r a = 0.0288/7 u c ReT 0 2 , and the thickness, 

S = OMzRc: 02 , becomes r w = 2.18Rc: 01 z/u t . 
Here the local Reynolds number is defined as 
R z. - Note that the coefficient c r is 

found to be a weak function of the local Reynolds 
number. Thus, for simplicity, we assume that c is 
constant. For Reynolds numbers corresponding to 
typical hybrid operation* c* is estimated to be 
approximately 0.55. In a real hybrid boundary layer 
with combustion and blowing, c- can be different 
from this estimation and for that reason, it is 
determined empirically. 

For the purpose of this paper, it is 
convenient to consider an average boundary layer 
delay and replace the local distance z, with die half 
of the grain length, LI 2, Note that the empirical 
constant c T accounts for die correction to the 
inaccurate selection of the length scale, L , 
However, we recognize that in reality there is a 
range of boundary layer delay times that should be 
considered. Thus the significant observation here is 
that a relatively broad band of oscillation 
frequencies is expected as opposed to a very sharp 
peak at the center frequency corresponding to 
Li 2. 
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accounted for by simply inserting time delays in 
the heat-flux expressions derived under ihe 
assumption of quasi-steady response. The 
implementation of this idea in the linearized 
version of Eq, 3 yields 

2l (*) = °*2 G\ (J ~ r bl \ ) — CTji?i (? - T bf2 ) (5) 

where 

t I T a » r Wl = T bJ\/ T t! > *b!2- T Mz/ T tl 

and r f/ « x/r 2 . 

Here r bn and r W2 arc the time delays 

experienced by the wall heat flux ( Q c ) to the 
changes in the oxidizer mass flux and the 
regression rate, respectively. The scaling of the 
time delays T bn and r^ /2 obey the general scaling 

law given by Eq. 4. However the c coefficients 
forti M and T bl2 are expected to be different since 

each of these delays represents a different 
adjustment mechanism for the boundary layer. The 
model presented by Equation 5 is central to the 
analysis in that it relates the heat conduction in the 
solid to the boundary layer combustion. 

6) Thermal-Combustion (TC) Coupled System 

Now, with the use of Equation 5 the 
thermal lags in the solid can be coupled to the 
combustion transients in the boundary layer. This 
coupling yields the following transfer function 
between the oxidizer mass flux (/(j) ) and the 

regression rate and represent the Hybrid 

Combustion block in the overall transient model 
(Figure 2), 


In our studies the response of the 
boundary layer to the changes in the mass flux is 


r TC 


u *>(*)_ 


2E Eg cr 2 e~ W1 j 


(A _ 'MW _ ^ Li-Zi— i- * 

f(s) (l + Vl+4j) Js+E^J-2E^ o +2E e j {e l +o- l e^' 2 j‘ 


( 6 ) 


This transfer function (TC coupled 
system) that represents the combustion 
phenomenon includes the dynamics of the thermal 


processes in the solid and approximates the 
combustion dynamics in the turbulent boundary 
layer of the rocket motor. Equation 6 can be used 
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lo investigate die stabfify c ha r acte r of the TC 
coupled system. The co nmoiily used method is to 
map the poles of the tnsfer fonction in the s- 
plane. Paikohriy, the real componen t of a ccrtaa 
pole of a transfer fvtdioa mdkales a mpfc fica tooa 
tMc associated with tot pde. Similarly the 
imagina ry part r e pres e nts the oscillation freq u en cy. 

First a system with zero bomdsry layer 
delays (Le. r M = =0 is de tenniae d to be 

always s ta bl e , la feet plot of the transfer fenrtan 
for a system wife typical pa r a m et e r s shows ao 
poles in fee s plane (only a zoo at (0, 0)) iadkatiag 
no sgn of instability (Le. Sec Figure 3). Ifi 
positive delay is ratrodracd betweea the ir.gir.witB 
rate and wall heat flax, a scries of fstabic poles is 
genera t e d. The example care with a 38 msec defay 
is showa in Figare 4. Note that afl the other 
pa rameters are kept i d ratirj l to the case wife ao 
delays. Even tho ugh there exists an mfinile scries 
of poles generated (at the same a mpl ifica tion rale) 
we wffl only concentra te on the pole with fee 
lowest frequency (fundamental mode). We believe 
that the higher order modes are an artifact of asfag 
a simple delay instead of the fed dynamics for the 
bo u n da ry layer traniarnK Mta w a ti the hi gh er 
frequency mo des, even if they exist m a read 
system, are fikriy to be Haipcd more effectively 
compared to fee fundamental mode. 

It is import rat to identify fee necessary 
conditions that mast ok for the ge n eratio n of TC 
coupled instabilities. A carefel namfea of fee 
denanxmafeor of the kaohr faction, Eq. 6, shows 
that for the instability to exist b l ow i n g ripe— I , 
4, ac ti va tio n energy, E m and delay tine, 
must be douzciol This mdirates that fee instability 
is a result of the coup li n g of three physical 
ph en ome na : Mocking of heat transfer fay radial 

solid and boundtey layer d ynamics The gas 
phasefeofid phase coupling mrrttaanm captured by 
the TC coupled model is shown in Fijpirc 1. 

Next we co ns i d er the effect of various 
parameters <m the TC coupled iasSabfety.lt can be 
shown feat the effect of the regi e l fera rale mass 
fhra e xpo n en t a, heal of vapor izati on of the fed, 
L r , the beat capacity of fee solid fed, C\ and fee 
surface t empe r alm c, T jy on fee osejOatioa 

fitqocacy and a mp B tad e we negligible for fee 
ran g e of these p aramet ers that are cocranonly found 
in hybrid systems. Also TC coupled system 
characteristics are completely independent of fee 


solid density, pj and the heat diffo&ivity in the 
solid, ar. 

The effect of the activation energy on the 
amplification and frequency is shown in Figure 5 
for 1*2=38 naffaeconds and k =0 j68. As 
mdirated by fee fipae, for activation energies 
commonly observed as hybrids, ranging from 5 
kcaVmole to 60 kcaVmole, fee vrafetion hi fee 
osci lla tion frequency is relatively small. The 
aupGIitation increases wife increasing adivajoQ 
energy. For this specific case the systems wife 
activation energies la rger than 3 kcal/mole show 
positive a mp lifkatk m aad u ns tabl e behavior. 
Decreased stability at ki^facr activatio n energies is 
ex p ected since at high E a values the s u r face 
temperature and regression rale are mere closely 
coupled. This coupling is ooe of the necessary 
mgir rfi rw K for ihr- pt nr rM ma nrf the instabilities. hi 

fed h the extreme case of E m =0, for which the 
smfece temperature is completely independent of 
the regression rate, oscillation frequency goes lo 
hxfhuty and the am p lificati on goes to negative 
infinity. 

The activation energies for the polymeric 
systems that are typically used in hybrids ran g es 
from 10 kca l ftnole to 60 kcafewle. For fee 
oonpolyTnenc paraffn-faused fuel system fee 
activraion ener gy is equivalent to the latent heat of 
vaponzteioa 1 *. This value for fee paraffin-based 
fuel formulations » estimated to be 17 kodfaoJe 
(Ref. 13). 

Similar results can be n hamrit for the 
effect of the blocking exponent, as shown in Figure 
6 for ^=038 and E. = 15 baVmolc. For the 
range of values repo t te d in the i t e rant for k (Le. 
0.6»4>.r7) lu *, the effect of blowing exjnmf oa 
fee frequency is neg E g iMe. The amplification 
increases wife increasing k and systems wife k 
values 0.45 and lager showed nestable behavior. 
Note that for fee physkaOy unrealizable erase of 
4=0 the ^stability d is a pp ea r s , since one of the 

finally Eke to note feat even though fee Mocking 
effect in paraffin based systems is somewhat 
reduced due to fee two phase character of the flow 
field, its is still significant enough to establish the 
coupling b etweea fee gas phase and solid phase. 

The most m flucn tiai p arame t er on the TC 
coupled frequency is detennmed to be the 
boundary delay tune, r^. Figure 7 shows the 
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predicted frequency as function of the delay time 
for three activation energy values, E m =5, 1 5 and 
50 kcal/mote. As shown m the %ae,lhe effect of 
the activation energy on the ft e q nrncy , espe cia lly 
an longer delay times, is smalL For the im p ose s of 
rim paper the effects of activation energy and 
Mowing opo n cgt wil he ignored and the 
fr equenc y wiD he represented as a function of the 
boundary delay time alone. The foflowiog carve ft 
rqnaticni is Tuggrrtrri to quantify the unrene 
relationship between rim frequency and the 



As shown in Figure 7, rins expressam 
represents the pre d i cte d ftag n acy to a very hqh 
degree of accuracy and will be used for all Ae 
hybrid systems tint wiD be discussed fader in tits 
paper. 

7) Cm Itynanrir Martrtr 

In the p re vi o us section tie dynamic 
models for the thermal lags and the hybrid 
combustion are mwarize d. In a d dition, rinse 
subsystem m o dels are c oupl ed to ob tain tie 

in the nans fln in the port. Ito ut m; firtir a fly 

as the c h a mber pressure; sp ec i fic imp u lse awl 
thrust cannot he o bta in ed solely from the TC 
coupled system. These variab le s can only he 
determined after introducing a model for the pt 
dynamics which will he discussed in tins section. A 
gas dynamic model, m general, should use tie 
oxidizer mass flow rate and tie fuel ns 
g enera ti on rate an the inputs aud it sho ul d yield tie 
parameters that are more clo s ely relat e d to tie 
per f orm a nce of rite motor such as tie cha mbe r 
p le asu r e and tic motor CVF ratio as die outputs. 

In car specific model (2V-pcrt model) we 
divide tie rfcnirr into tint componats* pis 
combustio n dumber, port volume wad tie posh 
comb ust i on chamhrr (Le. see Figpre K). We treat 
the pre and post combustion chamber s ns zero 
dim e as k na l voknne dements, but consider tie 
variation of state variables along the axis of tie 
port. The imthfinatic al ibnnobtiou is ob t ai ne d 
after the application of the cxmservatioo laws for 
the three volume el em e n ts . The set of equations 
resentifiug our model is solved analytically after 
liuearuation and as a result of that o p e r ati on, a 


transfer function for tic gas dynamic subsystem is 

obtained. Some m a n ei kal simulations are also 
pea feo ned to d et er m in e the validity of tie 
approximations used in the pertuhntion solutions. 
Note that a detailed description of tie gas dynamic 
model is given in tie Appendix of this paper. 

The following conclusions can be drawn 
from the results of the gps dynamic model as given 
in detail in tie Appendix of tins paper and in more 
detail in Rc£ 13. 

Gas dynamic model ca p lines tie 
filtig/en ytyhg dyna mi cs and also the 
lo tgilMdi nai a cou s ti c response of dm hybrid 
combus ti on chaatiri. As can he seen fr o m 
Figure 13 ac ou sti c inodes are stab le (well 
da mpe d ). No unst a b l e oscillation modes are 
generated by rise gas dynamic system by itself. 
The oocijncir du tabm et ootc the chamber 
acoustic modes. ‘Tiara! nil an showed tit a 
sudden pulse in the wiitforr rams flow rate 
caches the aco us tics modes. These modes 
decay in a relatively short period. Tbs ringing 
p h enom enon shows tint rim lo ngitudin al 
acoustic modes can be driven by nonfmear low 
frequency oscillations. Namely the nonlinear 
low fr equency pre s s are wa v e form s 
continuously cache tie acoustic modes of tie 
chamber. AMROC motor DM-Ol showed M 
least 5 kragitratm al a coustics modes. AD those 
were weft dam p ed . The gas d ynam i c model 
when corrected for iscatxopic speed of sound 
p red i ct s the a cro stics fr e q uen ci es of the DM- 
01 motor to a high level of accmcy (See 
Table 7). 

Gas dynamic modd coup l ed with the 
b o un dary layer dynamics (r M *0) does not 
produce any instabilities (no extra poles are 
generated in tie s pfane). This conclusion is 
based on transfer fcnrtioa analysis aud also 

hi g hlig h ts rim import rocc of the ti rimal lag 
system (Le. and its unique dynamic capability 
of prod u ci n g a phase lead) in p rodu cing a 
conp i rag mcrhmisni to gen e ra t e km 
frequency unstable modes. 

S ) Ther^Cun tin itl an G nu dyraunk (TCG) 

In this section we integrate the gas 
dynamic system to the TC co uple d system to 
establish the TCG cou pl ed response (Le. See 
Figure 2)l The TCG coupled system yields the 
per fo r m ance of the rocket (such as chamb er 
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pressure) for some given input of oxidizer mass 
flow rale, which is a fundamental control 
parameter m hybrid rockets. It is fair to state that 
the TOG coupled system represents the most 
fundamental dynamic behavior of a hybrid rochet 

The schematic of the hybrid subsystems 
and the a ft* na t i on flow between the sabsystems 
in the context of TOG coupling we shown is 
Figure 2. Note dot the lower block in the 
s ch em Mic represe n ts the TO coupled system. The 
input for the TO co u p l ed system is he local mass 
flu information and the o u tp rt is the mass 
ge n era ti on or the r egre ssi on rate. The u pper block 
shows the 2V-poit gas dy namic modeL The gas 
dynamic model takes the oudizer flow rate and the 
mass generation rate information and reveals the 
important per forma nce parameters such as the 
chamber pressure and specific impulse. The TOG 
coupled system can be considered as an overall 
model for the dynamics of a hybrid propulsion 
system dm utilizes gaseous oo o d iz e r winch is 
delivered by an isolated feed system. 

The generalization of the TOG coupled 
system to a liquid hybrid with vR i ifrat feed 
system d ynamics am be a chieved easily. Na mely 
two more modules for the Squid droplet 
ev a po rati on and die feed system dynamics be 
added to the front end of the TOG model. We 

systems does not involve either of Ae evaporation 
lags or the feed system dy nam i cs . 

In the TOG coupled system the gas 
dynamic module converts the fuel mass osciBaPons 
produced by the TO coupled system ito Ike 
chamber pressure osciRarious. It has been 
determined that an the p ro cess of conversion, the 
oscillations produc ed by the TO co u p l ed system 
(ix. fre q uenci es and a mp Cficati ou rates) are not 
altered. This fact is de m o nstra ted in Figure 9, 
winch shows a plot of the TOG system transfer 
function for the TO coupled system depicted h 
Figure 4. The same pk* also jmfcrUrs poles 
associated with the fflmg/emptymg mode (fiat 
order ncn-oscihiory mode) and he kapudad 
acoustic modes with negative real components 
(damped). TOG coupled system capites al the 
necessary tmdcriymg dynamics (fnemized version) 
of a gaseous hybrid with a de coup l ed feed system. 


9) NASA Ames Paraffin-Based Motor Tests: 

A promising fast burning paraffin-based 

fuel has been tested extensively in the Hybrid 
Combustion Facility (HCF) at NASA Ames 3 . The 
motor oabide diameter is 10 inches and gaseous 
oxygen is the oxidizer for all the tests c onducted at 
the HCF. The pressure time trace of a typical 
pwaffm-based hybrid motor test is shows in Figne 
10, whereas the Fomier transform for ihe same test 
is gives in Fignre II. Note dot the chamber 
p re ss ure traces for classi c al hybrids as pres e i n ed in 
the liter at ure are very samtbr m nature. As 
indicated mFignre 10 fora particular case, the feed 
system and the rhwb cr were com plet el y 
decoupled for al paraffin-based motor tests used in 
this paper. 

The Fomkr V ans form shows three broad 
peaks c orre sp onding to die low freq uen cy hybrid 
instability, the Hel mholtz mode and the first 
longitudinal mode. It has been generally observed 
that the do minam mode, low frequency osdBatioos 
are acco mpanied by lower amplitude higher 
fre q u ency acoustic modes. See Ref 6 for an 
extensive ( tisc u ssio u of the stability characteristics 
of the para ffi n based motors. 

If) Ci^trirw af Theury Results with the 
Mater Test Hutu 

In this section we wfll c o m p are the TO 
coupled theory mrgtahra frequency predictions 
with the chamb e r pressare osciftteioa fr eq u e n c ies 
from hybrid motor test data. In order to p roduce a 
conclusive c o m pari s on, we have selected a fange 
num ber of tests form several hybrid rocket 
development prog r ams . The key features of each 
p rogra m is discmsrd briefly in the following 
paragraphs. 

AMROC Motor Tests 443 : In this case the 
reported results are for 4 d ifferent thrust class 
hybrid motors which al ntilized LOXAGOX as 
the oxidizer and H i rtf as the fbeL The data 
from the AMROC tests is given in Table 2. 

■ HFDP 11 Ik* Motor Tests 5 *: The ckm&zct 
used in these tests is LOX md GOX and the 
fud is HTFB/Escorez focmabtiou The motor 
case d ia meter was 1 1 iaches and the d*a from 
5 tests used in ten sandy is given hi Table 3. 
Only the tests that exhibit c omb us ti on 
instability and p erforme d with single port 
motors are considered. The runs with the 
rearward-facing steps at the port entrance are 
also excluded m this study, since m these tests 
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the flow field in the port is significantly altered 

w hi c h may ca»c a unjor cfaagc m he 
boundary layer delay time coe ff icien t . We abo 
nclode test #1 of II inch LOX motor data in 
the analysis. Wc only consider test 1 since his 
die only ran with a chamber oonfi^ndn 
idetokal to GOX motors used in he 

- JIRAD Motor Tcste u:n : Ik oxidizer toed k 
these tests is GOX ad tbe fbei s 
HTHUEscorcz formbhoL The motor case 
ifia mctcr was 11 inches and the daeu firm 8 
tests mod in this study is given m TaUe 1 
Only the tests that exhibit c ombustion 
instability and per fo rmed with single port 
motors are considered here. 

Arizona Stale Umcfhy^ One test b 
m e f aded in he comparison for a 5 inch 
diame t er ni t ro n s oxide motor. The fuel was 
HTPB (ijc. See Table 5)- 
PXraffm-Based Molar Tests 9 : These tests art 
conducted a NASA Ames in a 10 inch 
dnmntikct test faeffity. The oxidizer was GOX 
and the foci is two separate nonpolyincfic 
para f fin -based formnfamons (Sf-li and Sf-4). 
The data form 25 motor tests are rep orte d in 
Table 6. 

The p re di c tion of the freq u enc y req ui res 
the csdnaiioB of the boundary layer delay time for 
each test from Ac rep o rted data. Since Ac reported 
data b so m e w ha t differen t for each p rog ra m 
d iffere n t scaling laws for the delay time mast be 
derived. 


The first one of those is for the various 
size AMROC motors winch operated at different 
L levels. Since the hdbnnadoa on the AMROC 
motors that cau be found an the ihciafturc b Imrilcd 
to L y c* and motor O/F ratio, it is desirable to 
express the boundary layer defay time, Eqnatioo 4 
(with z~L/ 2) in terms of those variab l es . The 
average velocity ni the port can be apprm in tor das 


With the use of the total mass flow relation, 
m t = and the defimriou L* = , 

the dday equation can be written as 

“ v m J?r ** ' 


Here V F /V m k die nfio of die port 

volume to the motor volume; which is (Jimrtid to 
be a p proi, wito cly OX for AMROC DM411 motor. 
We assume fot art AMROC a n to w s possess the 
same average gas oiiftwl temperature prodoct, 
RT^ and volume rtofo, - After noting dial 

all motors operate to very similar OflF ratios and 
values, it can be staled that n the series of 
tests the bouadnry layer characteristic dday time is 
proportional to I*. 

hi coder to reduce data form the rest of the 
tests, a similar relation for the bouadary layer dday 
time in tains of the rhaiuhrr p ressure and port flux 
levels must be derived. The port velocity can be 
r eplace d by an a v e ng e value (G„ + G t )fep mr 
with p^ given by the gas law P € = p m JRT mr 
to yield finally for the hciiwdnry layer lag: 


ur 


l bi2 


(C.+G,)W^ 


(ID 


Equation 11 gives the scaling law for the 
boundary layer defay m terms of the o pera tin g 
points of the rocket motor and also whh the size of 
the motor. This stales that dday in c reas e s with 
MC ft nag chamb er pressure, port length and 
d ec re a ses with mncpiig mass flux in the port. An 
important conclusion that can be drawn from those 
ob serv a ti ons b that m inc rea se in the cham ber 
p re ssure or a d e c re ase its the port mass flux act to 
decrease the rwcilfation frequency of this type of 
TC coupled system. 


„ g,Ki+ 2 o/rW+o/Fftcr,, 

~ 2 P c 

where is an average value in the port 

Substituting this expression is the delay fbranfa 
and using the idations for the total mass flow nte 
m t = G,A p and port volume V p = LA p yields 


The total mass flax can be rlaniatorrt form 
Equation 11 by introdneing flic average oxidizer to 
fbel ratio of the motor. 





V _A 

r,fo + 2G/FV(l4 O/FjjltT^ ' 


W Equation 12 wiH be used to reduce die paraffin- 
based motor data. 


9 

Ameri c an Institute of Aeronautics and A strotiauU cs 


A1AA-2003-4463 


Note that equations 10, 11 and 12 are 
equivalent mad die coBtatf mokipfier c is 
identical for all of these apicsskns. 

The avenge temperature gas ocnsM 
product dm a pp ea r s ia the draominrtor of afl foe 
bonadav layer delay tint conations is *— ^il to 
be comm for all of foe LOX/GOX motor tests 
that is consi d er ed ia dm paper- Tin is a faariy 
good assumption since KT is a weak f rac t io n of 
the motor Of¥ for the pra c tical opening 
conditions 0 . In fact the mniimi expected 
variation on foe RT valacs for all the LOX/GOX 
tests used for c o gy aiw n is predicted to be le» 
than 5% The ahsohne valne of RT m is selected 
to be 639 10 s (mfcecf far LOX/GOX motor tests 
and since the same valne is esed for aH calcafarin 
any error in RT^ wifl only effect the anooerical 
value of the empirical delay constant c . 

In order to doable check the KT ^ vahe 
used in the orientations one could use the fast 
acoustic frequ e n cy air aw e d form the motor test 
data and nse tet to rslana t r Ike speed of soaad 
and finally KT^ insed on the following equation. 


Here / is the ratio of the specific heats averaged 
over die length of the motor, fy_ L is the first 
lougituririiad acoustic mode and L, is the effective 
length of the motor. For a move acc urat e 
calculati on the cc pph i gas dynamic model 
discussed in die Appendix of thk paper can be 
used. Ad of the NASA Ames motor tests with 45 
inch long gai ns had hear first acoustic mode at 
around 370 Hz. For these tests with nse of y*125 

and £*=1.27 m, one ob tai n s KT^ of 7.1 10 s 
(m/sec) 2 which is only 10% hig j h e t than the 
B urned value KT^ for the Up test is se l ected 
to be 5.1 10 5 (m/sec) 2 . This is 20% l o an dam die 
RT^ for the GOX/UOX system due to die low 
te m per a t — c of foe N^Odiydrocarhoa combustion 
pro duct s . 

Now Equation 7 can be coupled with the 
expre s s i ons for the boundary delays time to 
evaluate the frequencies predicted by the TC 
coupled modeL 


f\~L ~ 


JFkt. r 


<13) 


/ = 0.48- 


KT^ 


Ki>o/fV(i+2o/^ 7Z, 


(14a) 


0,48 

clP r 


I ' 

G m KT^ 

Oih\ 

cLP € 


(14b) 


Tfie oscaflation fieqnency as a fu n c tion of 
the boundary delay time predated by die TC 
c ou ple d theory a l ong with data form foe hybrid 
motor development programs are plotted is Figure 
12. Note font we have nsed a hwnday layer delay 
cons ta n t value, c , of 2J01 for afi motor tests h 
order to obtain the best fit between foe theory and 
test results, ftenc note foal fins t tprrimrn ta l vahe 
for c is the same order of magnitude with foe 
rough est im a t e fix. —0.55) for a boundary layer 
with no blowing and no combustion. 

Considering the wide range of test 
conditions, the agreement b e twe e n the measured 
pressure oscillation frequencies and foe TC 


cou ple d model predictions is surprisingly good. 
Nole foal the comparison mrfartrs 43 motor tests 

• Using resnks from 4 differe n t 
devtlo pmcm and resea r ch p r ograms 

• Using 3 separate oxidizers (Le. LOX, 
GOX andN^Q) 

• Covering a wide range of motor 
dim en sio n s (Lc. from 5~ OB Arizona state 
motor op to HT OD AMROC motor) 

• Using many d iff ere nt fuel compositions 
(Le. HTFBs, HTTWESCOREZ, 
HTPB/DDI/ESCOREZ, 2 pmffia-hnscd 
fonnulathxis) 
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• Covering a wide range of operating 

conditions. 

Equation 14a indicates that for motcn 
operating at the similar O/F, the oscillations 
frequency is inversely proportional to the L of die 
motor. Tins conclusion is conskteat with the 
previous results reported in the btodnt. Simibriy 
Equation 1 4b shows for motor ope r atin g at equal 
oxidizer mass flaxes, rhmhtf jicssaes and OF 
the oscillation frequency dc c i cne s with aacvensng 
fuel gram length , ta order to prove dial motor 
length plays a critical role setting the asedtuboa 
frequency, we have plotted the Ames motor data, 
wirida con ta ins two seb of grain lengths, by 
renxwing the length effect hum the scaling law. 
This resulted in m mcrease in the scatter in the 
fr eq uency data indicating the impo rta n ce of the 
grain Wwgth- A sanrihr ch ec k to cwfi r ni be 
importance of length can be conducted by 
c o m p arin g the HRAD and HFDP minor tests (Lc. 
ran with grains 102 inches long) with the Ames 
paaffe tests (Le. ran with g rams 32 or 45 inches 
long) that were conducted under very spnabr 
oper a t i n g c o n riitio ux, For example comparing test 8 
of the HTDP progra m to pmafii u famed motor test 
4P-01 shows dnt gram hag* uant he i n cl ud ed n 

in ihc oscOhtion fr eq uencies for these tats. The 
same con clu sio n can be <kawn by comparing test 2 
of HFDP program to the test Ames test 4L-08. 

TC coupled theory predicts oscittatious of 
the re gre ssi on rat c/ fuc l mass generation rate aad 
would be observed directly an a hybrid motor as be 
radial osci lla ti ons of the d iff u si on frame within the 
boundary layer. The gas dynamic model converts 
the fuel mass genoteon rate oscillations to the 
chamber pressure oscillations. It is expected that a 
system operating at higher OF ratios wadd 
p ro du ce less rha a nbtr pre ssu re oscillations sace 
fuel mass c oustihrte a smaller fa ction of the told 
mass expelled from the nozzle. This fat » 
observed in NASA Ames motor tests, namely togh 
CVF motors, in general, presented lower mrplitades 
c hamb er puw r t oscillations cmya red to the 

The sample linear model p r es ented in this 

amplitude low frequenc y oscilfaions and be 
scaling of the frequency. But it abo falsely predicts 
mi indefarte growth of cha mber pressure 
oscillations. In reality, the nonlinear mechanisms 
that exist in a hybrid motor would limit the 
indefinite growth of the csdltatious and result in 


limit cycle oscillations that are commonly observed 
in motor tests. Moreover the simple theory does not 
explain when these low frequency instabilities will 
take place and what their amplitudes would be. h 
has been shown in the hybrid development 
p rogr ams that the fare end co ufign rub m/volmne of 
the motor pfays * paramount impo rta nc e in setting 
die amp lit u de of the low frequency instabilities. 
The motor tests indicated that, at least for gaseous 
motors, axial injectors result m more stable 
operte k m compared to the radial o conical 
mjectiou of the oxidizer. 

A natural e xplana tion for these 
observteions is dnt the TC coupled instability 
m e chanism almost always exist in hybrid systems, 
bat the amplitude of the Emit cycle oscillations will 
depend on the excitation level by the disturbances 
m the right frequency range. We bcheve that the 
fore end fluid dynamics (Le. vortex s h eddi ng ) is 
die pru na ry source of di st u rb a nces in a hybrid 
rocket. For exanqdc changing die injection scheme 
of the oxidizer from mad to radial, may introduce 
a flow distu rb a n c e component dtat k preferred by 
the TC coupled system. If one assumes th* all the 
low frequency oscillations are developed by the 
fluid medumics at the fore end of the molar, one 
would fa! to expfan the length effect on the 
oscillation frequency. Also afl oxidizer np ccti on 

seem to produce the same oscillation freq u ency. 
This fat is abo efiffkull to explain by die 
assumptions that afl the esedfarions are produced 
from the fluid dynaunc events at the fore end of the 
motor, since the injection scheme should after the 
fluid n rech a nirs s ignificant ly (Le. vortex shedding 
frequency). 

We would fike to note that most of 
musical instruments also work on the same 
principle 2 *. For example in a flue organ pipe the 
column of ar inside be pipe is set into vibra tio n 
(Le. at its natural f req u encies ) by an edge tone 
which b produced by a jet of air that i m ping es on a 
carefully de si gned Ep. The vortex sheddag M the 
lip sets the edge tone which itself excites the 
natural inodes in die pipe. Note dm in be case of 
an organ pipe the jet speed and be Bp geometry is 
carefully selected for the pipe lengh such that the 
edge tone produced at the fip ooatrius a large 
component at the nafinl frequency of the pipe that 
needs to be cached. In this organ pipe malog, the 
bp is the fore end of be motor (Le. geometry, 
injection sche m e, velocity etc) and the organ pipe 
is die hybrid system m od el ed by the TC coupled 
theory. 
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11) Conduskms: 

The following conclusions can be drawn 
from the results of the TCG coupled model: 

A phasible mec hanism that geacnlcs km- 
frequency diamber pressure oscillations is 
fflg g rUrd . This mechanism b due to the 
coupling of the thermal lags, the gas phase 
combvtioa and gas dynamic sub-systems of 
the hybrid rocket The physical transient 
node! of the motor is presented in a 
ma th ema t ical formalism that allows one to 
rtfm al r Ac e xp ecte d oscOlatioB frequency. 
The amplitude of die oscillations can not be 
determined by fob linear modeL 
The oscillation frequency estimated by the 
model b in very good agreement with the 
motor test data from several p r ogram s utilizing 
three dis tin c t oxidizers* several fuel 
fonnutations, a wide range of motor 


growth of the oscillations. In reality the 
amplitude of foe chamber pressure fluctuations 
will be limited by nonlinear effects that are not 
covered in the linearized TC coupled theory. 

It has been shown that the longitudinal 
acoustic oscillations that are commonly 
coexist with the low frequency oscillations can 
be generated by the ex ci ta tio n induced by foe 
nonlinear waveforms of the low frequency 
c hamb e i disturbances, in short, it b plausible 
that high frequency acoustic modes are driven 
by the low frequency osculations. 
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scaling law for foe primn y hybrid oscHbtim 
fr eq u ency: 

/= 0239 (2+— (15) 
l O/F ) LP C 

639 10 ?{mhecf far GOX/LOX system 
XT mv = 5.11 10* (mfaec) 2 for low energy oxidizm 
systems such as NtO 

This mode of testability b common to all khd 
of hybrids: liquefying or co nv enti onal fuels, 
liquid oxidizer or gaseous CKjdzzcrs. f! b 
reasonable to believe that fob low frequency 
mode b present m every hybrid system to 
some rrtert Seme motors Me m ore unstable 
compared to foe others hujaa these coolers 
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of the chamber pressure, it faih to establish an 
estimation process for foe ampfifcnfe of the 
oscillations. The theory predicts an unfimited 
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14) Appendix 

In our gas dynamic model the motor is 
divided nto four subsections: 1) pre-combustion 
chamber, 2) port volume, 3) post-combustion 
chamb er, 4) nozzle. Each component is 
approximated app ropriately with one or zero 
spatial dimensional models. The subscripts 1 and 2 
are used for the physical properties in the pre and 
post comb u stion chamb er s r e sp e cti vely. Various 
stations of significance are also distinguished by 
use of subscript notation. Namely the entrance 
station of the gaseous oxidizer into the system is 
represeAed by a subscript o, the nozzle entrance is 
demoted by n and the port entrance and exit are 
shown by a and h, re spec ti vely. Our models for 
these subsections are discussed in the next section. 

Pre-Combustion Chamber 

We use a zero dimensional gas dynamic 
model for die precombustion cha mber . Thus the 
momentum equation and the energy equations 
reduce to the mean stale i nfo rma ti on of the 
pressure and temperature in this volume. But, the 
mass balance equation is nontrivial and it can be 
written as 



Note that the size of the vohme is 
assumed to be uaitat over the hn tune. We 

polytropic process. Peep*' in the cha mber. Here 
n f is the potytropac exponent We b e liev e that 
these are reasonable a ssu m pti ons for ibis relatively 

combus ti on- Under these stated comBtiom and with 
the use of for ideal gas bw, cquotiou A1 cn be 
reduced to 

dP y _ P, ^ 

aj = «*„ m a wheir o,= (A.2) 


13 

Am eri ca n Institute of Aerorarics and Astramrics 


AIAA-2003-4463 


Port Chamber 

In our model we treat the port as a quasi- 
1D tube wkh continuous mass and heal addition 
The conservation of mass can be written as the 
following difi a qta l equation- 


dp d G 

dt dz J 


(A 3 ) 


Note that for the simplicity of notation we treat the 
gas co ns tant tt nye ntm t product as a single 
dependent variable. 

Furthermore we assume that the gas 
constant temperature product at the exit of the port 
volume is a fraction of the nozzle entrance value 
RT e which caa be ra l rnlrtrd from the equilibrium 
chemistry at the o p e ra! M g O/F ratio of the motor. 


Here mj is Ac mass generated at the axial 
location z of the port. At das stage we prefer to 
express die mass g en e r atio n ten ■ Ac general 
f u nc ti onal form. 

»/ =Ag* z ) 

The m o ment um equation in the port caa be written 
in terms of the local mass fhu, G „ the local 
pressure, P aad the local density, p . 

££. fkiA LL { (; / c K (A „ 

d t d z d z y 2 A p J p 

Note fal the last term ■ thu apatuni a 
included to fjpfttc he effect of dtc skin fr i ctio n oa 
the force balance. Here Cy stands lor fre skin 

friction coefficient a location z, A r as fre poet 

area and the C for the ckc u nfeitact of the port. 
All the geometrical p roperties are assumed to be 
independent of time and axial dimension. This ba 
reasonable a ppro xima tion for typical hybrids which 
posses a uniform port diameter ifatointkn afrng 
the kngdi and also a slow variation of the port 
geometry 15 . 

We further a ss u me dm the ideal gas bw 
holds locally in the port. 

p=/*r (AS) 

The closure of die system requ i res tike 
energy equ atio n. For the sake of siuyfki ty we 
replace die energy equation with the following 
linear variation of temperature gas constart pnrkt 
in the port. 

R7\z)= (KT^ + - (*?•), W* <**> 


*T 2 =/ k RT c (A.7) 

hi short, we specify the whole temperature 
field in die motor chamber for a gives operating 
point and we i gnore the variations in the 
temperature field dung the t i au rim l s to that 
respect 2 V-port is an c o t bumal gas dynamic 

— - j-t 

mood. 

Past-Combmtkm Chamber 

Similar to the pre-combustion chamber 
ir ramMri , we assign uniform therm odynamic 
p rope rti e s for the gas in die post-condmstion 
cha mb er volume. Th erefore the dynamics can be 
repre se nted with the mass balance which takes die 
following form after the fixed volume assumption. 



Here m h is the mass flow rale entering the hot 
volume aad m m is the mass flow rale exiting the 
volume through the nozzle. Note that we have 
inherently assumed that post-combustion chamber 
behaves like a well stirred reactor with the 
reactions fast compared to the transient rakes of 
concent. Equation A_8 can be written as 

dP V 

a 2 -^-=m k -m m where a i=-jfiT- <*•*) 

We have a ppr oxi maled the gas constant 
t empera ture product in die post-combustion 
cham b er with produc t at die port exit. 

Nozzle 

We first assume that the nozzle flow is 
quasi-steady which is only valid at very slow 
transients. We also assume that duongh oat the 
operation of the motor, die nozzle is always 
choked, namely at any tune the nrio of the 
chamber pre ssu re to the ambient pressure is larger 
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than the critical value associated with the 
pro p e rti es of the expending gas. After these 
smpfifyhg assumptions the nozzle flow Me can 
be expressed with use of the cl ass i c a l quasi- ID gas 
dynamic formula 






+*. 


_d 

dr 


(All) 


Here k m represents the effect of inertia of the gps 

m die c o nvergen t portion of the nozzle in caging a 
phase dif ference bet we en die oscOMon of the 
velocity aid pressure at the ertiance of the nozzle. 


where 



Here y € is the ratio of the sp e cific heats of the 
combustion products. 

At moderate frequen c y oscillations (sodb 
as the aco usti c modes of the ntoi^ the acoustic 
a dmitta n ce of he nozzle most be rowiadnt d and a 
more s o p hi s tica ted behavior must be m od el ed. For 
our purposes, it is adrqnatr to ada p t Crocco's 
expression 34 for the nozzle transfer function which 
is valid for small dev iati on s from die mem 
chamber pr essur e and for mod erate frequenci es * 


(A.H) 

V 2 Cm) 

This integration is per forme d on the 
whole couvogent portion of the nozzle. Note that 
Cm is the critical sound velocity, a is the local gas 
velocity in the nozzle, v r is the gas velocity at the 
nozzle entrance. 

Equations A .9 aid All and the ideal gas 
law can be combined to relate die p re s sure 
pert u rbation directly to the mass flow rale 
p erturba ti on al the exit of the port. 




(A 13) 


Note that a fraction of the nozzle volume 
is added to the post-combustion chamber volume to 
account for the inertial effects of the gas in ibe 
nozzle. For typical hybrids the effective nozzle 
volume » small co mp are d to the post-combustion 
chamber volume. We finally indicane dot for oar 
isothermal tr eatme n t die temperature derivative 
term is identically zero. 

Equation AI3 hrhrlct some properties 
of die co m b u stio n product s w hich are not known 
offhand. These p rope rti e s ae calottatrrf from 
chemical equi lib ri u m conside rati ons with use of 
the code STANJAN for the selected propeQatts. A 
sample calculation result of the product, RTc with 
the motor Of ratio for HTPB-cnygen system is 
given in Ref. 13. 

to this section we develop snail 
perturbation solutions and also nwnwiral solutions 


for the system of expiations derived for the 2V-port 
model . 

Perturbation Solutions 

We duen mn t perturbation solutions 
anxmd a n ominal op e ratin g pond of the rocket for 
the gas dynamic system of equations described ra 
die fvevious paapqhs These perturbation 
solutions wifl be nsed to obtain a transfer function 

sy s t em . The traisfer function nses the o tifag 
mass flow vale as die input and die system p re s s ur e 
for the output The details of the derivation are 
shown in Ref. 13. 

Tbe transfer fraction specifically b et we e n 
the oxidizer mass flux input and the pre- 
combustion chamber p re ssu re cutfwl can be written 
as 
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0> " (U-, +^) *+*] A+ktF*. -«.M x+(2L,fr^* (l 2 ,*e)|l«.k(A£/2) ' ' 


Tins function alsoicpvesesbftc ovtnfl m t ac of 
the response for gs dynamic systen. Sumhr 
expressions far Ae other variables can be also 
derived, bat far the impose of ov arguments An 
expression k a d equa t e. 

^V w m r^ 5 hato}iu 

hx outer to ch eck the validity of the 
assumptions used in the perturbation solutions, we 
also perform laanrricaf mbtioas on tbc 2V-port 
model described in the previews section. We stui 
with the fink differencing of the special d airaiw 
terms and bsc tbc second order central difluuaj n 
over a mxfbnn mesh along Ac axis of the port 25 . 
For the tone m a r cha a g we nse the RwpsKsti 
fourth order method. The boundary oo o dh ow for 
the port equation s arc ob ta ined firms the pre and 
post co m b us ti on dumber dynamic cqnarinm The 
numerical method used in the c a i r n ialir u s gave 
stable and accurate results lor a typical that step if 
3x10* seconds and 80 mesh poin ts in the port. The 
details of the mncrical scheme is discuss ed in Ref 
13. 


The simulation code is applied to various 
transient events such as the throttling and the 
stnosoidal oxidizer flow rate case. Some of the 
results are di s cus s ed in the following section. 


The respon se of the gas d y nam i c system 
for a hybrid motor can be deduced from the p& 


stability features of Ac system arc revealed by the 


which is Ae rharretcriitk cqnnkw However my 
in v e st iga ti on of Ae gas dynamic system r e spo ns e 


ic tf o a ne s an input for the cornbesfiop re sp o nse Am 
is represented by the term. 

in this section we nse the simplest quasi - 
steady co w iws tkwi modd. Mnwrfy , dl Ac tr an sien t 
effects inchitfmg Ae Aermal fags in Ae sofid and 
also the bomadary layer delays except Ae gas 
d y namics are ig n ored h this modd. For Am 
reason, this modd which does not dorr any 
dynamic coupling b et w ee n Ae sub-systems of the 
rocket is highly de si r ab le to isolale Ae 
investigations to purely gas dynamic phenomena. 
The combustio n respon se f un ctio n cm easOy be 
d eri v ed front Ae g en end recession rale law after 
lin e ariza t i on and Laplace trandormatioa. 

F n =a,»G*~ , (A.15) 

For Ae case of Ae quasi-steady 
combustion Ae gas dy na mi c system t r ans fe r 
function is plotted m Fignre 13. The operational 
paramrlns of this motor are dnenssrd in Ref 13. 
The plot shows a pole on Ae negative real axis and 
a row of imagwary p o le s also h Ae negative real 
inn pone, ivc real pofc icscjiidkj me Hwuncuit 
capacity of the chamber, wh e reas Ae imag i nary 
poles re pre s ent the tongitnrfinal acoustic nature of 
the cavity. Note that all the poles have negative 
read parts, namely the boA acoustic and also Ae 
global response components decay in time. The 
c o n cl u sio n is that Ae pure gas dynamic system is 
stable m nature. 


In order to test Ae validity of Ae transfer 
function we consider Ae trivial case of no aims 
addition and no l onpe rat u arc or velocity gradients 
io Ae port. A tins case Ae transfer function 
simplifies to 


B x = fim 


JRT^J 


*h) -.M 

V 7 V ~ a +) 


(A.16) 
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After the application of the L’HospftaTs rule one 
obtains the well known transfer function for a 
simple single volume system. 

e =_ Old- (a.17) 

(^+12+1) J + e TjJ + l 

Hoe Tj = L fe* is tie characteristic filhng time 
of the motor chamber. The other prancters arc 
defmed as £ =V[A, and c which 

are era— only used m rocket h tm Ct. Note tint 
the total v ol—ir of the moior co mbc ai oa cha mb er 
can be written as f = (^+^+1^ .Thsiith 

transfer function generally used (e sp e cia lly in the 
industry) in hybrid dynamic modeli ng 


The p re c e ding equation does not predict 
the filling/emptymg time scale of a hybrid rocket 
accurately. For that reason we develop a better 
estima t e by the fo B oa ag argument, hi the region 
of the s plane that the pole for the ffllmg/emptying 
process b expected to be locded, the following 
apixuumatioiis can be made. 


Aa K-M (^i*) 

Under these approximations the gas dynamic 
tran s fe r Sanction anp ifi c s to 


{t, +L 7 +ji , -Lj -fl^jL/ 2 ) M+e ^ - unh^/^ 

Note flat since - a ^ term is independent of s, this transfer function repres ent the behavior of a first 

order system with a characteristic time of 


Li + 



M 

T f ~ 

e 

jl-ttsfcpVr-a 

m 



<A20) 


The throttling response for a tfaruttlhg 
ntio of 12 of the motor with the w 
spe cificati ons is ca lc ula ted with tfifltitm m et h ods 
and plotted in Figure 14. Figure shows the tine 
variation of the non-dimeasiooal pre-corabastion 
c ham ber p re s sure defined as 

pressur e and P/ is die final p re ss u re. The striking 

observation is dkat the perturb a tion sohboa with 
the fifing tame defined by tqaaiwi A-20 is in 

resnlL This sh o w s thn equation A20 gives a good 
r^iimar for the fllhg the of the motor c ha mber 
within the context of the 2V-pcei model. The 
re sp onse acco rdin g to the simple transfer f un c ti on 
r% also shown m Figure 14. h is obvious from the 
figure that die simple filling time alim en t for the 
hybrid combustion chamber is not valid. 


Another observation from Figure 14 is 
that the nu meric a l sanitation result shows 
oscillatory be havi o r in the early stages of the 
thrott l in g event. Has osdftatory behavior b due to 
the excitation of the acoustic modes of the cavity 
by the disturban ce caused by the sadden increase in 
the oxidizer mass flow rate. The emanation 
frequency and the riaayiag time scale of this mode 
of motion is coauWal with the longifnitinrl 
acrosti c mode fan lings for the same motor. 

This kmd of be ha vior is co mm o n ly 
observed in AMROC DM-Ol motors 4 . Nandy, a 
nonlinear disturb— cc causes the wcD darned 
oscillation of the c hu nb er pte saa e at the 
kmgjhnfifial aco ustic frequ e nci es . The anting 
disturbances in the act— I tests are bel ieved to be 
caused by: 1) a sadden bl ockage of the nozzle 
throat area by some dg—egraied fad fractions, 2) 
the ignition of a bulk of fuel oxidizer mixtur e 
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captured in the pre-combustion chamber or 3) the 
low frequency oscillations of cham ber pressure. 

We c o m pa re the predictions of our gas 
dynamic model for die acoustic frequencies of the 
chandler with the AMROC DM-O! motor test data. 
Table 7 shows the first five longitudinal 
frequencies d e t er min ed by the pertnrtaxn 
solutions. Since our model is isothermal, these 
frequencies are corrected by multiplying with tie 

factor of Jy . An avmgr y vdv of 12 is used 

for die calculations. Note that the c orrec te d 
frequencies falls in the range of observed 
frequencies for every cilmhirri aconstic mode. 

The gas dy nam ic code dtitkyed to 
perform simula t i ons or the 2VTut model is used 
to confirm the findi n gs of fir putuilutieu 
solutions on the frequency respo n se of the system. 
In the simulalionsy sinusoidal oikfcncr mass flour 
rale inputs at various frequencies are applied as the 
forcing function and the phase and the anphtude of 
the dumb er pressure oscillations are determined. 
The agreement b et ween Ac prilmlntiun sofcxtiocs 
and the numerical si— ilatiom for the faquwy 
re sp ons e of die system is also found to he 
satisfactory. 

The simulation code is alio used to 
investigate the response of the gas dynamic system 
to a ( fetutim cc. The ouidire r mass flow input is 

1 5) TaMes ami Figures: 


suddenly increased 10% above its operating value 
and hdd at its elevated level for 0.(13 seconds 
before its suddenly reduced to its original value. 
The pressure time history for this simulation is 
shown Figure 15. It is clear that the disturbance 
excites the acoustic inodes of the cham ber similar 
to the previous observations. These exited modes 
decay in a moderate time scale as it was observed 
in the AMROC tests. A similar simulation is 
perfonned to determine die effect of nozzle 
blockage on the stability by red u c in g the nozzle 
t hro at area for a brief period of thne. The results 
for the nozzle area perturbation case me i d raikal 
in nature to the resnfcs of the oxidizer mass flow 
rale perturbation. In short, a nonlinear disturbance 
canes the ringing of Ac hybrid cham ber which 
decays in time. 

Finally a simple delay (r M ) on the 
c ombus ti on r esp o nse is applied to the gas d ynam i c 
system. 

Fjr = a^nG^'e'™ (AJU) 

Both die p erturba ti on solutions and also 
the numerical smnbtiom indicate thU any positive 
combustion delay have negligible effect on the 
stability nature of the system. 





1 ) Solid phase kinetic times 

T v < 10' 3 sec 

Degradation mrrhum nn of the 

2) Gas phase kinetic times 

r_<10' 3 see 

HyAucaxbon conduction 

3) Feed system response 
times 

(Varies greally from system to 

Response time of the feed system 

4) Evaporation tiroes 


Evaporation process of the liquid 
axabaaor 

5) Thermal lags in soSd 

r* ocr/r 2 aHT 1 sec 

Thermal profile changes in fie solid 
#aiu 

6) Boundary layer diffesion times 

%«!/»,* 10 1 see 
(Varies greatly form case to case) 

Tnrhnirni boandary layer difinrion 

3) Aconstic times (longitudinal) 

r a at Lie ® 10 -3 see 
(Varies isedy form case to case) 

Propagation of the acoustic waves 

7) Gas dymunc tiffing times 

Tjgg ac L/c * JO 1 sec 
(Varies greatly form case to case) 

Global mass flow balance 


Table I: The transient time scales of various phenomena in a typical hybrid rocket. Only the last four arc 
modeled in this paper. 
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ESS3 

OtF 

• 

(ft/scc) 

• 

f a«p 

(ft/sec) 

L* 

(“) 

/ 

(Hz) 

T m 

(msec) 

S Motor 
(10k) 

1.45 

5584 

5361 

570 

11 

48.9 

Half- Scale 

03k) 

1.45 

5584 

5361 

1305-1430 

4.8 

117.4 

H-500 

(75k) 

1.45 

5584 

5361 

1770-2340 

4.0 

1763 

DM-01 

(250k) 

155 

5700 

5472 

2168 

2-3.5 

204.8 


Table 2 : Summary of paramr i r js used in the frequency estima ti ons for AMROC motors. 


Test 

No- 

D 

(m) 

(■*ec) 

G. 

(Win 2 

«) 

G, 

(fc/ in’ 
«) 

K 

(p«*> 

O/F 

/ 

(Hz) 

A^ e 

(!»*) 

r «2 

(msec) 

1 

(GQX) 

4-25 

6.0 

0-423 

0571 


2.85 

13 


488 

2 

(GOX) 

4.25 

m 

0.402 

0541 

550 

289 

8 


47.2 

7 

(GOX) 

420 

m 

0339 

0.458 


289 

15 

100 

44.7 

8 

(GOX) 

4,22 

5.1 

0365 

0.492 

260/280 

287 

20 

H 

255 

i 

(LOX) 

* 

• 

0.222 

031 

435 


H 

250 

66.2 


TaUe 3: Su mm a ry of parameters used is the frequency es tim ations for the 1 1 inch Hybrid motor tests* 


Grain length is 102 inches* 


Test 

Nol 

D 

(«) 

L 

(in) 

■. 

(Wsec) 

G. 

(IMn 2 - 

sec) 

G, 

(wW 

B-SCC) 

K 

0*0 

OIF 

/ 

(Hz) 

<%) 

T U1 

msec) 

3 

Ksa 

102 

3.44 

038 

0527 

745 

Ub 

6-10 

15 

665 

4 

43 

102 

336 

023 

0325 

325 

2.4 

10-20 

20 

47.4 

6 

3.7 

102 

836 

0.78 

1.024 

750 

32 

10-15 

12 

33.7 

7 

53 

102 

3.41 

0.15 

022 

335 

21 

6-15 

15 

73 3 

8 

El 

102 

387 

015 

mm 

mm 

2Jb 

2-5 

60 

191.7 

9 

6.1 

102 

634 

028 

0392 

215 

2Jb 

620 

33 

25.9 

15 

H 

108 

1020 

0.45 

0599 

1025 

El 

4 


83.9 

10 

EH 

102 

5.93 

020 

0278 

215 

2.6 

10-25 


365 


Table 4: Swnnary of parameters used in the frequency es tim a ti on s for the NASA/MSFC 1 1 inch Hybrid 
motor tests* 


19 

American tastitule of Aeronautics and Astronautics 








































































































































AIAA-20Q3-4463 


Test 

N* 

D 

<-> 

L 

<«) 

(ftvscc) 

G. 

(fc%r- 

sre) 

(ps) 

0/F 

/ 

(Hz) 

T W2 

(msec) 

1 

2.0 

27 

3.44 

039 

£JjJI 

4jD 

19 

20.6 


Table 5: Summary of parameters used in the frequency es timati ons for the Arizona Stale Hybrid motor 
lest 


Test Nil 

L 

0") 

G. 

(Arin^-sec) 

Pc 

(P®*) 

OIF 

/ 

(Hz) 

T U 2 
[msec) 

4F-4 

32 

0.44 

528.0 

3.97 

41j6 

133 

4F-5 

32 

0.49 

551.0 

339 

39.9 

136 

^lb 

32 

0.20 

561.0 

372 

14.4 

29.7 

4F1-C 

32 

0.16 

5430 

3.06 

13j6 

36.1 

4F-01 

45 

039 

318.0 

369 

40.8 

123 

4P32 

45 

038 

9933 

348 

133 

383 

4P-03 

45 

0.40 

939.1 

365 

14.4 

35.6 

4L-03 

45 

031 

641.6 

369 

137 

30.9 

4L-04 

45 

032 

656.7 

366 

23 .6 

18.9 

4L-05 

45 

0.46 

6493 

372 

19.1 

213 

4L-08 

45 

0.44 

525.0 

2.64 

233 

17.9 

4MH 

45 

038 

318.7 

2.40 

39.6 

134 

4P-04 

45 

021 

159.1 

1.73 

434 

10.7 

4Lr09 

45 

036 

2653 

134 

203 

13.7 

4L-10 

45 

0.43 

590.0 

389 

273 

183 

4L*1 1 

45 

0.11 

2133 

136 

17.8 

273 

4L-12 

45 

0.13 

3013 

301 

1534 

332 

4NF-01 

45 

039 

6023 

377 

2360 

23.1 

4NF-02 

45 

032 

6004 

334 

11.10 

404 

4NF-03 

45 

036 

500.8 

396 

2330 

21.4 

4NF-04 

45 

0.48 

5683 

331 

3375 

183 

4L44 

45 

0.43 

5243 

351 

25.15 

183 

4ST-02 

45 

0.45 

5403 

371 

2436 

183 

4L-15 

45 

031 

5553 

2.09 

14.94 

263 

4Rep-02 

45 

0.45 

4023 

363 

36.49 

133 


Table be Sammy of parameters used in the frequency estimations for the NASA Ames Hybrid motor 
tests. 
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Eataated 

CNRdri 

Observed 

FltiMSlld 

AemntieMode 


freyty 

Firfirry 

AmpfifiestM 


(Htr) 

(Hz) 

(Hz) 

(1M 

i* 

48.4 

53.0 

50-55 

-22 

r* 

96.7 

103.9 

100-110 

-22 

y* 

142 j 6 

1 156.1 

150-165 

-21 

4* 

191.1 

2093 

200-220 

-19 

5 * 

240.7 

263.6 

250-275 

-18 


Table 7: Thr esti mate d aid observed longjtodinal acoustic frp qacBQ O for AMROC PM-01 motor. 




Figwt 2. Tbc schematic of the TCG coupled system. 
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Figure 5: Efl^ erf aefrvabcn eneigy on the c Bcflbt ioB freqaeacy and anpfifiotioBL 



Figirt ft: Effect of blocking factor expmcnt on the oscillation frequency and amprifkatknL 
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Figure 7: Effect of boundary layer delay tiine on the oscillation frequency for various activation energies. 
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Figare 9 z The contour plot of die TCG coupled system with = 38 msec. Most of the transient features 

of a hybrid system can be deduced from the locations of the poles. 



Figare It: Chamber and feed pressure tine traces for the paraffin-based motor test, 4K)5_ 
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Figure 11; Fourier transform of the chamber pressure for 



Fipro 12r Comparison of the TCG coupled oscillation 
data available in the liter a t ur e. 
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H gre 15c The Tcspomse of the prp-combctkMi chamber prrssmr to a dcaulxm e ofthccBodcer 
flow ole. 
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